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ABSTRACT 

We use a suite of cosmological, hydro dynamical simulations to investigate the chemical 
enrichment history of the Universe. Specifically, we trace the origin of the metals back 
in time to investigate when various gas phases were enriched and by what halo masses. 
We find that the age of the metals decreases strongly with the density of the gas in 
which they end up. At least half of the metals that reside in the diffuse intergalactic 
medium (IGM) at redshift zero (two) were ejected from galaxies above redshift two 
(three). The mass of the haloes that last contained the metals increases rapidly with 
the gas density. More than half of the mass in intergalactic metals was ejected by 
haloes with total masses less than 10 11 M and stellar masses less than 10 9 M . The 
range of halo masses that contributes to the enrichment is wider for the hotter part of 
the IGM. By combining the 'when' and 'by what' aspects of the enrichment history, 
we show that metals residing in lower density gas were typically ejected earlier and by 
lower mass haloes. 
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1 INTRODUCTION 

Despite the fact that 'metals' - heavy elements produced 
by stars - make up a very tiny portion of the cosmic matter 
budget, they are of critical importance to our understanding 
of the Universe. From a diagnostic point of view, they are 
very useful in constraining star formation and tracing feed- 
back from massive stars and supernovae. More importantly, 
however, metals impact the rate at which gas cools, affect- 
ing structure formation on scales of galaxy clusters down 
to dust grains. This effect should not go understated since, 
depending on the temperature, metals can reduce the cool- 
ing tim e by over an order of magnitude for solar metallic- 
itv (e.g.lSutherland fe Dopitall 19931 : iGnat fe Sternberdl2007l : 
I Wiersma et al.ll2009al ~ 

Observations clearly show that a substantial fraction 
of the diffuse intergalactic mediu m (IGM) has b een en- 



l2010l ; IPieri et al.l 12009). However, the physical mechanism, 
the timing and the sources of the enrichment all remain 
unknown. The IGM may have been enriched by the first 
generations of stars and galaxies at very high redshift s, 
or it may have been polluted by more massive galaxies 
at intermediate redshift s. Metals may be carried out into 
intergalactic space by galactic winds driven by supernovae, 
radiation pressure from star light, or by active galactic n uclei 



(AGN) (e.g ICen fc Ostrikei] Il999 



i Madau et al.1 l200ll: iTheuns et al 



Aguirre et all l2001allbl: 
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20041: lAguirre et al] 120051; lOppenheimer fe Pavel 12 006: 
Dave &; Q ppenhei merl 120071: lOpp enheimer & Dave 
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riched with heavy elements (e.g , 



Schav e et al.l 
Simcoe et al 



Elliso n et al.1 
Aracil et al.l 



Cowie & Songaila 



2000 
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Schav e et al 
Schave et al 
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2007: 



Wiersma et al.1 l2009bl: IShen et al.1 120101 : iTornatore et aD 
2010l ; IChoi &; Nagaminell2010l ). In dense environments such 
as groups and clusters of galaxies metals may also be 
efficiently dispersed via processes such as ram pre ssure 
stripping or tidal stripping (see e.g. Bor gani et al.ll2008l for 



Aguirre et al.1 



Danf orth fe Shulll 120081 : ICooksev et aD 



a review). 
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The most straightforward way to investigate when the 
IGM received its metals, is to plot the metallicity or metal 
mass as a function of redshift, as we did for a series of cos- 
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mological simulations in I Wiersma et al.l (|2009bl ) . By com- 
paring models with different physical implementations, one 
can obtain information about the processes that are impor- 
tant for IGM enrichment. However, the metals that reside in 
the IGM at high redshift may not be the same metals that 
are in the IGM at the present day. What is hot IGM gas at 
high redshift, may have cooled by low redshift. Intergalactic 
metals may also accrete onto galaxies and end up in the ISM 
or in stars. 

Here we also use cosmological simulations to study the 
enrichment of the cosmos - noting that the IGM is of par- 
ticular interest - but we take a novel approach in that we 
classify the gas based on its state at the redshift of inter- 
est (we choose z — and z — 2) and trace its enrichment 
history back in time with the aim of understanding how the 
observed metals were put into place. Specifically, we will ad- 
dress the following two questions: When did gas that resides 
in the IGM at z — and z — 2 receive its metals? And 
what are the masses of the haloes that last contained these 
metals? 

This paper is organised as follows. In Section [2] we de- 
scribe the simulations used, emphasising the details that are 
particularly pertinent here. Section [3] discusses when the gas 
was enriched, taking each sub-phase in turn. In Section 3] we 
investigate what enriched the gas, using a halo- finder to de- 
termine what kinds of objects are mainly responsible for the 
enrichment of various gas phases. We connect the two ideas 
in Section [5] relating enrichment time with the last halo that 
contained the metals that end up in the IGM. Finally, we 
summarise our findings in Section [6] 



2 SIMULATIONS 

The simulations presented here are part of the Overwhelm- 
ingly Large Simulations (OWLS) project (| Schave et al.l 
2010), which consists of a large suite of cosmological, hy- 
drodynamical simulations that include chemodynamics. We 
use a substantially modified version of the iV-body Tree- 
PM, smoothed particle hydrodynamics (SPH) code gad- 
get hi, which is itself a modified version of the code gad- 
get II ([Springell l2005h . The code includes new module s 
for radiative cool ing and heating ([ Wiersma et al.l [2009a]) , 
star for mation (| Schave &; Dalla Ve cchial l2008h , chemody- 
namics JWiersjna^^iJ_|2009b|), and feedback from super- 
novae (JDalla Vecchia &; Schavd 2008). Below we will give 
a brief overview of the simulation ingredients that are 
most relevan t for t his work. More details can be found in 
I Schave et al.l (|2010l ) and the references given above. 

We assume a flat ACDM universe with the 
cosmological parameters [Q m , Q& , Qa , as , n s , /i, Y ] — 
[0.238,0.0418,0.762,0.74,0.951,0.73,0.248], as determined 
from the WMAP 3 -year data and co nsistent^ with the 
WMAP 7-year data (|jarosik et al.ll2010h . 

We use t h e rad iative cooling tables presented in 
I Wiersma et al.l (|2009al rl. In brief, net radiative cooling 
rates are computed element-by-element in the presence 



1 The most significant discrepancy is in as, which is 8% lower 
than the value favoured by the WMAP 7-year data. 

2 We used their equation (3) rather than (4) and CLOUDY version 
05.07 rather than 07.02. 



of the cosmic m i crowa ve background (CMB) and the 
lHaardt fe Madaul (J200lh model for the UV/X-ray back- 
ground radiation from quasars and galaxies. The contribu- 
tions of the eleven elements hydrogen, helium, carbon, nitro- 
gen, oxygen, neon, magnesium, silicon, sulphur, calcium, and 
iron are interpolated as a function of density, temperature, 
and redshift from tables that have been pre-computed using 
the publicly avail able photo -ionization package cloudy, last 
described by Ferl and et al.l (|1998l ), assuming the gas to be 
optically thin and in (photo-) ionization equilibrium. 

Stars are formed from gas particles as per 
ISchave fe Dalla Vecchial (|2008h . Gas at sufficiently high 
densities (tih > lQ ~ 1 cm~ 3 ) is expected to be multiphase 
and star- forming (|Schavdl2004r ). Because we lack the reso- 
lution and the physics to model the cold interstellar phase, 
we impose an effective equation of state with pressure 
P oc pj for densities exceeding tih = 0.1 cm -3 , normalised 
to P/k = 1.08 x 10 3 cm -3 K at the threshold. For this slope 
of the equation of state both the Jeans mass and the ratio 
of the Jeans length and the SPH kernel are independent 
of the density, thus preventing spurious fragmentation 
due to a lack of numerical resolution. Star- forming gas 
particles are stochastically converted into star particles at 
a pressure-dependent rate that is determined analytically 
from the observed Kennicutt-Schmidt star formation law. 

Supe rnova feedback is injected in k inetic form, as de- 
scribed in iDalla Vecchia &; Schavd (J2008) . After a short de- 
lay of £sn = 3 x 10 7 yr, corresponding to the maximum 
lifetime of stars that end their lives as core-collapse super- 
novae, newly- formed star particles inject kinetic energy into 
their surroundings by kicking a fraction of their neighbours 
in a random direction. The referen ce simulations presented 
here use the default parameters of IDalla Vecchia &; Schavd 
(2008), which means that each SPH neighbour i of a newly- 
formed star particle j has a probability of r]rrij / ^2-^ h rra 
of receiving a kick with a velocity v w , where Af ng b = 48 
is the number of SPH neighbours. We choose r\ — 2 and 
v w = 600 kms -1 (i.e. if all baryonic particles had equal 
mass, each newly formed star particle would kick, on av- 
erage, two of its neighbours). Assuming that each star with 
initial mass in the range 6—100 Mq injects 10 51 erg of kinetic 
energy, these parameters imply that the total wind energy 
accounts for 40 per cent of the available kinetic energy for 
a Chabrier IMF and a stellar mass range 0.1 — 100 Mq (if 
we consider only stars in the mass range 8—100 Mq for 
type II supernovae, this works out to be 60 per cent). The 
value r\ — 2 was chosen to roughly reproduce the pe ak in 
the cosmic star formation rate fsee lSchave et al.ll2010l ). 

For t he chemo-dynamics, w e employ the method de - 
scribed in I Wiersma et all (|2009bn . We use a lChabrier] (|20Q3h 
stellar initial mass function, with mass limits of 0.1 Mq and 
100 Mq. We consider enrichment due to mass loss from AGB 
and massive stars as well as type la and II supernovae (SNIa 
and SNII respectively). Our prescription for stellar evolution 
is heavil y based on the Padova models as we use the life- 
times of IPortinari et al.l (119981) and t he yield s of low mass 
and h igh mass stars of iMarigol ([200 ll) and IPortinari et ahl 
(|l998h . respectively. The yields of IPortinari et al.l (|l998h in- 
clude ejecta from SNII along with estimates of mass loss 
from high mass stars. Using these yields gives an element of 
consistency between the low a nd high mass star s . Fin ally, 
we use the 'W7' SNIa yields of iThielemann et al.l (J2003T ). 
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As the precise channels of SNIa progenitors are uncer- 
tain, we employ SNIa explosion rates that roughly reproduce 
the observed cosmic SNIa rate. We note that the observa- 
tions from different surveys are not consistent and that our 
rates conform to the most recent data. An e-folding delay 
time is used to describe the SNIa rate and we use a corre- 
sponding efficiency (i.e. number of white dwarfs that become 
SNIa) of 2.5%. 

We use SPH-smoothed metallicities for cooling and stel- 
lar evolution purposes, where Z sm = pz/pg for gas parti- 
cles, and star particles inherit the sm o othed metallicities 
of their progenitors. In IWiersma et al.l (2009b) we argued 
that this definition is more consistent with the SPH for- 
malism than the more commonly used particle metallicity, 
Zp a rt = mz/rrig. Using smoothed metallicities results in a 
spreading of metals over greater volumes, increasing the av- 
erage rates of cooling and star formation. 

Since SPH is a Lagrangian scheme, metals move 
with the particles they are attached to, and hence a 
high-metallicity particle may in principle be completely 
surrounded by low- metallicity gas. Introducing smoothed 
metallicities partially mitigates this numerical artefact. 
However, SPH codes may underestimate the true amount 
of ga s mixing as compa red to grid cod es in case of instabil i- 
ties (JAgertz et al.ll2007h or turbulence ([Mitchell et al.ll2QQ9h . 
It would thus be very interesting to investigate whether 
the metal distribution in such grid simulations is similar 
to t hat obtaine d with SPH. On the other hand, we note 
that IS hen et al.l (|2Q10h included turbulent diffusion of met- 
als in cosmological SPH simulations but found its effect to 
be small. 

Of the OWL S suite of simulations (described in 
ISchave et al.ll2010l h we use the subset listed in Table [1] Here 
we make use of the notation LXXXNYYY, where the sim- 
ulation cube has a comoving side length of XXX /i _1 Mpc 
and uses YYY 3 dark matter particles and YYY 3 baryonic 
particles. Table [1] also lists the dark matter and the ini- 
tial baryonic particle masses, radm and rab, the maximum 
proper gravitational softening (the softening is fixed to this 
value in comoving coordinates above z = 2.91), e pr op, the 
redshift at which the simulation was stopped, z en d, and a 
short description. 

We explore two variations on the reference model 
(REF). Run AGN_L1 00N512 includes the prescription of 
iBooth &; Schavd ((2009) for the growth of supermassive black 
holes by gas accretion and mergers, and for feedback from 
AGN. The black holes inject 1.5 per cent of the accreted rest 
mass energy as thermal energy in the surrounding interstel- 
lar medium. This model re produces not only the the local 
black hole scaling relations (Boo th &; Schavd boog). but also 
the observed entropy, temperature, and metallicity profiles 
of groups of galaxies an d the stellar masses that they contain 
([McCarthy et alJl2009h . 

For run WDENS.L100N512 the parameters of the 
winds are dependent on the local density, keeping the energy 
per unit stellar mass formed constant. They scale such that 
in high-density environments the wind velocity is higher: 




-5 -4 -3 -2 

Log 10 (Mass weighted average of Z) 



Figure 1. Mass weighted metal distribution in temperature- 
density space at z = in simulation L050N512. The colour scale 
gives the metallicity. The contours indicate the metal mass distri- 
bution and are logarithmically spaced by 1 dex. The straight lines 
indicate the adopted division of the gas into: non-star-forming gas 
(i.e., n H < 0.1 cm -3 ), diffuse photo-ionised IGM (p < 10 2 (p), 
T < 10 45 K), cold halo gas (p > 10 2 (p), T < 10 45 K), WHIM 
(10 45 K < T < 10 7 K), and ICM (T > 10 7 K). The metals are 
distributed over a wide range of densities and temperatures. 



7] 



7?* 



n H 



-1/3 



(2) 



n H 



1/6 



(1) 



where ^ = 600 kms -1 and 77* = 2 such that the model 
agrees with the reference model at the threshold for star 
formation. For higher gas pressures the velocity increases in 
proportion with the local sound speed. ISchave et al.l (|2010l ) 
find that this model is much more efficient at suppressing 
star formation in massive galaxies, though not as efficient as 
the model that includes AGN. 



3 WHEN WAS THE GAS ENRICHED? 

In this section we will investigate how the time at which a 
gas element was enriched varies with its density and temper- 
ature at some later time. We will evaluate the enrichment 
history of the gas as a function of its physical state at both 
z = and z = 2. Before doing so, it is, however, instructive 
to look at the distribution of metal mass at the present time, 
which we show in Fig. [T] 

Fig. [T] indicates how we can divide the temperature- 
density plane in several regions of interest. For nu > 
0.1 cm -3 (which corresponds to p/{p) > 5.3 x 10 5 at z — 
and p/{p) > 2.0 x 10 4 at z — 2) we impose an effective 
equation of state and assume the gas to be star-forming. 
We divide the non- star- forming gas into hot gas typically 
found in haloes of large groups or clusters (ICM, T > 10 7 K), 
warm-hot intergalactic gas (WHIM, 10 45 K < T < 10 7 K), 
diffuse photo-ionised IGM (p < 10 2 (p), T < 10 45 K), and 
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Table 1. List of simulations. From left-to-right the column show: simulation identifier, dark matter particle mass, baryonic particle mass, 
maximum physical gravitational softening (the softening is held fixed in comoving coordinates above z = 2.91), final redshift, and a brief 
description. 



Simulation 


m dm 


ra b 




Cprop 


^end 


Comment 




(^Mq) 


(^M ) 


(h 


-^pc) 






L025N512 


6.3 x 10 6 


1.4 x 10 6 




0.5 


1.45 


Reference model 


L050N512 


5.1 x 10 7 


1.1 x 10 7 




1.0 





Reference model 


L025N256 


5.1 x 10 7 


1.1 x 10 7 




1.0 


2 


Reference model 


L050N256 


4.1 x 10 8 


8.7 x 10 7 




2.0 





Reference model 


L100N512 


4.1 x 10 8 


8.7 x 10 7 




2.0 





Reference model 


AGN.L100N512 


4.1 x 10 8 


8.7 x 10 7 




2.0 





Includes AGN 


WDENS.L100N512 


4.1 x 10 8 


8.7 x 10 7 




2.0 





Wind vel. prop, to local sound speed (SN energy as REF) 



cold halo gas (p > 10 2 (p), T < 10 45 K). While we have di- 
vided the cold (T < 10 4 ' 5 K) gas into low- and high-density 
components (diffuse IGM and cold halo gas, respectively), 
we have not done the same for hotter gas. We find that if 
we split the WHIM at p — 10 2 (p), the enrichment history 
of the low-density component is very similar to that of the 
WHIM in total, whereas that of the high-density component 
is more similar to the enrichment history of the ICM. The 
high-density component contains less than one quarter of 
the WHIM metals at z — and less than one sixteenth at 
z = 2. 

In I Wiersma et al] (|2009bh we showed that the metallic- 
ities predicted for these phaseqj by our reference model are 
in agreement with the observations, although the metallicity 
for the diffuse IGM may be slightly too low and is not yet 
fully converged with respect to the numerical resolution (it 
increases with increasing resolution). 

The contours in Fig. [1] show how the metal mass is 
distributed in temperature-density space. The metals are 
spread over a wide range of phases. The metallicity (that 
is, the ratio of metal mass to total mass), indicated by the 
colour scale, ranges significantly between the various regions. 
That the contours do not always trace the regions with the 
highest metallicities highlights the fact that the metals are 
distributed rather differently from the rest of the baryons. 
Table [2] summarises the z = and z — 2 metal mass distri- 
butions for the reference model. 

We will make frequent use of what we will term "metal- 
mass- weighted mean redshifts" , zz ■ This quantity, which we 
compute on the fly for all gas particles in our simulations, is 
defined as the mean redshift of enrichment events weighted 
by the metal mass imparted by each event: 



Table 2. Fractional distribution of metals in the reference models 
at z = and z = 2. 



zz 



J2i Amz-iZi 



raz 



(3) 



where Amz-,i is the metal mass received by the gas particle 
at redshift z\ (i indicates a time-step) and mz is the total 
metal mass of the particle. Thus, if a gas particle received all 
its metals over a short interval (which is typical for particles 
that have been ejected from haloes), this quantity gives the 
redshift at which it was enriched. Note that zz does not 
correspond to the time at which the metals were created. 



3 Note that in I Wiersma et al.l d2009bh we used T = 10 5 K as the 
dividing temperature for WHIM and cold gas; we lower the limit 
here so that we can better isolate the components visible in Fig. [2] 



Phase Fraction of metals (%) 


L050N512 (z = 0) 


Stars 


68.2 


SF Gas 


5.8 


ICM 


1.3 


WHIM 


19.0 


Cold Halo Gas 


4.6 


Diffuse IGM 


1.0 


L025N512 (z = 2) 


Stars 


42.0 


SF Gas 


25.8 


ICM 


< 0.1 


WHIM 


24.0 


Cold Halo Gas 


7.2 


Diffuse IGM 


0.9 



Stars may create metals long before they are released and 
metals may be cycled through multiple generations of stars. 
We will nevertheless sometimes refer to the metals residing 
in gas with higher values of zz as being "older" . In JH we 
will show that for intergalactic metals zz is typically close 
to the redshift at which the metals last resided in a halo. 

Fig. [2] shows zz in temperature-density space for 
L050N512 at z = (left panel) and for L025N512 at z = 2 
(right panel). Note that this figure contains no information 
about how much metal mass is actually in a given bin (for in- 
stance, the diffuse IGM contains relatively little metal mass, 
but this is not obvious from this plot). We can, however, 
refer back to Fig. [1] to see the metal distribution. For refer- 
ence, computing the metal mass weighted average over all 
gas (stars) yields 2 Z »1.5 (z z ~ 1.2) for L050N512 at z = 
and z z « 3.1 (z z » 3.2) for L025N512 at z = 2. 

The most prominent trend visible in Fig. [2] is a strong 
overall gradient with density, with higher density gas be- 
in g enriched more recent ly. This agrees with the predictions 
of lOppenheimer &; Pavel ((2009) for the low-redshift, diffuse 
IGM as traced by O VI absorbers. Of all components, the 
star- forming gas (i.e., the ISM) has been enriched most re- 
cently. This is not unexpected, as this is the gas that sur- 
rounds the stars that are releasing the metals. Moreover, 
much of the metal mass that was injected into the ISM a 
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Figure 2. Enrichment redshift in temperature-density space. The colour scale gives the metal mass weighted mean enrichment redshift 
as a function of the gas density and temperature at z = in L050N512 (left panel) and at z = 2 in L025N512 (right panel). Each 
pixel shows the metal mass weighted average value of zz for all particles that fall into the corresponding temperature-density bin. The 
rectangles in the bottom-right parts of the diagrams indicate the corresponding values for the metals locked up in stars. Also plotted 
are t coo \ = £h contours for solar (solid) and primordial (dashed) abundances. Note that we have stretched the colour scale to emphasise 
the different components (only a handful of pixels correspond to high values of zz) and that the two panels use different colour scales. 
Denser gas was typically enriched later. The enrichment redshift is particularly small for gas with a short cooling time. 



long time ago, will already have been locked up in stars. 
Indeed, the metals locked up in stars were, on average, re- 
leased (by other, older stars) earlier than the metals that 
reside in the ISM as can be seen by comparing the colour of 
the rectangles in the bottom-right parts of the panels to that 
of the ISM. The fact that the increase of zz with decreasing 
density continues all the way to the lowest densities, sug- 
gests that travel-time is a limiting f actor fo r the enrichment 
of the IGM, as suggested bv lAguirre et all (J2001bh . 

Fig. [2] also reveals some interesting trends with the gas 
temperature. For high gas densities, the dependence on tem- 
perature reflects a dependence on the radiative cooling time. 
We can see this by comparing to the solid and dashed con- 
tours, which show where the net radiative cooling timqj, 
tcooi, equals the Hubble time, £h, for gas with solar and 
primordial abundances, respectively. The two contours that 
enter the diagram at high and low temperatures correspond 
to net cooling and heating, respectively. They converge at 
high densities and intermediate temperatures at the equilib- 
rium temperature. There the contours should really extend 
further to the bottom right, approximately to the star forma- 
tion threshold (above which we impose an effective equation 
of state), but they stop due to the finite resolution of the 
grid on which we have computed the cooling time. Gas to 
the right of the contours has £ coo i < ^h , but gas to the left of 
the contours cannot cool. Clearly, the enrichment redshift is 
a decreasing function of the cooling time if £ coo i < £h- Since 

4 We plot only the radiative cooling time, but note that the ac- 
tual cooling rates in the simulation are dominated by adiabatic 
expansion for low densities. 



the gas cannot remain in the same area of the temperature- 
density plane for more than a cooling time, this implies that 
the time since enrichment is correlated with the time since 
the gas parcel moved significantly in the T — p diagram. This 
is for example expected if the gas was enriched and shock- 
heated at roughly the same time, as in a galactic wind. It 
would, however, also happen if the enrichment led to a strong 
reduction in the cooling time. 

In the long cooling time regime (i.e. to the left of the 
contours) the enrichment redshift also depends on the gas 
temperature. We will show later (see Fig. [3]) that some of the 
features for p < (p) and T > 10 4 K are sensitive to the imple- 
mentation of feedback. However, the finding that the coldest 
part of the diffuse IGM was enriched more recently than the 
WHIM, is robust. Note, however, that the diffuse IGM con- 
tains very few metals compared to the WHIM (see Fig. [TJ. 
We have checked and have found that most of the metals 
in the cold, underdense part of the diffuse IGM have never 
been part of a (resolved) halo. This suggests that this gas 
may have been enriched by winds from intermediate-mass, 
intergalactic stars (presumably formed in galaxies and then 
ejected). Apparently, gas that has been enriched and ejected 
by galactic winds ends up in the WHIM at higher temper- 
atures, probably because it has been put onto a higher adi- 
abat as a result of shock-heating. Once the gas has entered 
the region where £ CO oi > £h, it can no longer decrease its en- 
tropy, preventing low-density WHIM gas from moving into 
the diffuse, photo-ionised region of the diagram. 

The left panel of Fig. [3] shows zz(T, p) at z — for the 
L100N512 version of the reference model. Comparison with 
the left panel of Fig. [21 which showed the prediction for a 
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Log to p/<p> 



Log to p/<p> 



Figure 3. A comparison of the enrichment redshifts for simulations including different feedback processes. The colour scale gives the metal 
mass weighted mean enrichment redshift as a function of the gas density and temperature at z = for simulations REF-L100N512 (left 
panel), WDENS.L100N512 (middle panel), and AGN.L100N512 (right panel). Model WDENS uses an implementation of supernova- 
driven winds that is more efficient in massive haloes than that of model REF. Model AGN uses the same prescription for supernova 
feedback as model REF, but includes feedback from AGN. The rectangles in the bottom-right parts of the diagrams indicate the values 
of zz for the metals locked up in stars. For numerical reasons that are explained in the main text, the unresolved multiphase ISM (Log 
p/(p) >6) follows a much tighter equation of state for models WDENS and AGN than for model REF. While the features seen at 



(p,T)~ (10- 

are robust. 



(p) , 10 K) and (10 (p) , 10 K) for model REF depend on the implementation of galactic winds, the other trends 



simulation with 8 (2) times higher mass (spatial) resolution 
shows that the trends have converged numerically. 

The remaining panels of Fig. [3] show the enrichment 
redshift predicted by two other simulations from the OWLS 
project. The middle panel shows model WDENS, for which 
the wind velocity scales with the local sound speed, while 
keeping the total energy per unit stellar mass identical to 
that used in model REF. In this model the supernova driven 
winds remain efficient up to higher halo masses, which re- 
sults in a low er global star formation rate at low redshift 
(|Schave et al.l l2010h . The right panel shows the results for 
model AGN, which uses the same prescription for supernova 
feedback as model REF, but which also includes feedback 
from accreting supermassive black holes. Since AGN feed- 
back is most relevant for high halo masses, the star formation 
rate, and thus the rate at which metals are both produced 
and locked up in stars, decreas es much faster below z — 2 
than for the other models (Schaye et al. 2010). This explains 
the increase in the metal mass weighted mean enrichment 
redshift of the stars, from zz ~ 1.1 for model REF to 1.5 
for WDENS to 1.6 for AGN (the corresponding values for 
the gas are zz ~ 1.1, 1.3, and 1.5, respectively). 

The gas on the equation of state (Log p/(p) > 6), which 
represents the unresolved, multiphase ISM, displays much 
less scatter for WDENS and AGN than for REF. This re- 
flects the fact that for the latter model the output time 
corresponded to a global time-step whereas passive parti- 
cles were caught 'drifting' adiabatically between time-steps 
in the reference run, resulting in temperatures that deviate 
slightly from our imposed equation of state. Note that this 
does not represent any difference in how the hydrodynamical 
forces were calculated among the three simulations. Observe 
also that some of the particles in the right panel fall on a 
line near 10 8 K. This corresponds to the AGN heating tem- 
perature (see Booth &; Schavd 12009) and thus implies that 
these particles have been heated very recently. 

Comparing the three panels, we see that the features 
seen at (p, T) - (1Q" 1 (p) , 10 43 K) and (lO -05 (p) , 10 57 K) 



for model REF depend on the implementation of supernova- 
driven winds. All other trends of zz with T and p are, how- 
ever, robust. Even AGN feedback does not change the trends 
discussed above. Since the differences between the simula- 
tions are small compared to the similarities, we will only 
make use of our higher resolution reference model for the 
remainder of this section. It should be kept in mind, how- 
ever, that this model probably underestimates the typical 
enrichment redshift due to the absence of an efficient feed- 
back mechanism in massive galaxies. 

Figs. [2] and [3] do not allow us to infer what the typical 
enrichment time is of metals as a function of just density or 
just temperature, because the diagrams do not take into ac- 
count the distribution of metal mass in the T — p plane. The 
solid curves in Fig.|3]therefore show the metal mass weighted 
mean enrichment redshift as a function of density (left panel) 
or temperature (right panel) . The black and red curves show 
results for L050N512 at z = and for L025N512 at z = 2, 
respectively. 

Focusing first on the left panel, we see strong gradients 
with the gas density. While gas that is diffuse at z = was 
typically enriched at z ~ 2, metals that reside in collapsed 
structures (i.e. p> 10 2 (p)) at the present time, were on av- 
erage released much more recently. The enrichment redshift 
continues to decrease down to the highest densities. The 
same holds at z = 2, although the gradient is less steep at 
high densities. Gas that is diffuse at z = 2 typically received 
its metals at z > 3 and metals residing in underdense gas at 
z — 2 were on average ejected at z > 4. 

Similarly, the right panel of Fig.Ushows the metal mass 
weighted enrichment redshift as a function of the gas tem- 
perature at z = (black) and z = 2 (red) for all gas (solid) as 
well as for gas with density smaller (dot-dashed) or greater 
(dashed) than 10 2 (p). The enrichment redshift peaks for 
T rsj 10 5 — 10 6 K, but this merely reflects the more funda- 
mental trend with density: metals in this temperature range 
mostly reside in diffuse gas while lower and higher temper- 
atures are dominated by higher densities. 
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Figure 4. Metal mass weighted mean enrichment redshift as a function of the gas density (left panel) or the gas temperature (right 
panel) at z = for model L050N512 (black) and at z = 2 for L025N512 (red). The different line styles correspond to cuts in the 
temperature (left panel) or density (right panel) as indicated in the legends. The curves were cut when the number of particles per bin 
falls below 200. The enrichment redshift decreases rapidly with increasing density. Metals residing in gas with intermediate temperatures 
(T rsj 10 5 — 10 6 K) are typically older, but this merely reflects the fact that metals at these temperatures reside mostly in lower density 
gas than metals at lower or higher temperatures, as can be seen by comparing the different line styles in the right panel. 
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Figure 5. Cumulative probability distribution of the enrichment redshift, weighted by metal mass, of metals at z = (top row) and 
z = 2 (bottom row) for stars and for the different gas phases defined in Fig. [T] The different colours correspond to different numerical 
resolutions. The results are converged with respect to resolution except for the diffuse IGM and the WHIM. For this gas even our highest 
resolution runs may underestimate the enrichment redshift. There is substantial scatter in the enrichment redshift, particularly for the 
phases whose metals are, on average, older. 



So far we have only plotted the metal mass weighted 
enrichment redshift averaged over all gas particles in some 
temperature and/or density bin. Such plots do not reveal 
information about the distribution in the values of zz of 
the particles that fall into the bin. Fig. [5] provides this in- 
formation by plotting the metal mass weighted cumulative 
probability distribution function of zz- In other words, it 



shows the fraction of the metal mass that resides in gas par- 
ticles for which zz exceeds the value that is plotted along 
the x-axis. The different columns and line styles show the 
different regions in T — p space defined in Fig. [T]and the dif- 
ferent colours correspond to different resolutions. The top 
and bottom rows are for z = and 2, respectively. 

Fig.[5]shows that the metals residing in a particular gas 
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Figure 6. Metal mass weighted mean log 10 TTihalo: where ra^alo is the total mass of the halo that last contained the metals, as a function 
of the gas density and temperature at z = in L050N512 (left panel) and at z = 2 in L025N512. Each pixel shows the metal mass 
weighted average value of log 10 ra^alo f° r all particles that fall into the corresponding temperature-density bin. The rectangles in the 
bottom-right parts of the diagrams indicate the corresponding values for the metals locked up in stars. Only particles that have resided 
in haloes with more than 100 dark matter particles are included. Note that the two panels use different colour scales. Metals residing 
in lower density gas originate from lower mass haloes. Because high-density gas (p> 10 2 (p)) currently resides in haloes, the increase in 
halo mass with temperature reflects their virial temperatures. The diffuse IGM is nearly completely missing from this plot because it is 
enriched by unresolved haloes and/or intergalactic stars. 



phase typically span a wide range of enrichment redshifts. 
Note that this plot underestimates the true range, because 
it shows the distribution of zz for gas particles, which have 
themselves been averaged over all past enrichment events 
(weighted by metal mass). That is, even though a particle 
may have experiences multiple enrichment events, it only 
has a single value of zz- The distribution is particularly 
wide if the average enrichment redshift is high. While the 
average value of zz is about 2 for the diffuse IGM at z = 0, 
a quarter of the metal mass resides in particles with zz < 1 
and another quarter has zz > 3. 

Comparison of the different line colours shows that the 
convergence with numerical resolution is excellent for stars, 
the ISM, cold halo gas, and the ICM, but that it is less good 
for the diffuse IGM and the WHIM, particularly at z = 0. It 
is likely that our highest resolution simulation at this red- 
shift, L050N512, still significantly underestimates the typi- 
cal enrichment redshift of these low-density components. We 
will show in the next section that this is because low-mass 
haloes dominate the enrichment of low-density gas. 



determining the mass of the halo from which the metals 
were last ejected. 

We identify haloes using a friends-of-friends group 
finder with a linking length of 0.2. This group finder is 
run on the dark matter particles in the simulation. Bary- 
onic particles are associated with their nearest dark matter 
neighbour. We include only haloes with at least 100 dark 
matter particles, because we found that th i s limi t results in 
excellent agreement with the ISheth et al] (|200lh fit to the 
dark matter halo mass function. For model L100N512 this 
lower limit corresponds to a typical z = stellar mass of 
~ 10 90 M© and this is also roughly the value below which 
the stellar mass function cuts off due to the limited reso- 
lution. For models L050N512 and L025N512, however, the 
stellar mass already cuts off below about 10 2-5 dark matter 
particles, which corresponds to stellar masses of ~ 10 8 ' 3 M© 
for L050N512 at z = and - 10 77 M© for L025N512. We 
therefore expect that we will underestimate the importance 
of the lowest mass haloes that we consider, particularly for 
these last two models. 



4 HOW MASSIVE WERE THE OBJECTS 
THAT ENRICHED THE GAS? 

Now that we have a general impression of when the cosmic 
gas was enriched, we would like to investigate what types 
of galaxies were responsible for the enrichment. That is, do 
galaxies residing in certain halo masses preferentially enrich 
a typical phase of gas? We will address this question by 



To make the connection between metals and the haloes 
out of which they were ejected, we perform the following 
calculation. For each enriched particle, we begin at z — 
(z = 2 for the L025N512 simulation) and search for that 
particle (using its unique ID) in the friends-of-friends halo 
catalogue. If said particle is found, we note the total mass of 
the halo (rahaio), and the corresponding redshift, Zhaio- For 
all particles not found in this output, we choose halo outputs 
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Figure 7. As Fig.[6]but for simulations REF_L100N512 (left panel), WDENS.L100N512 (middle panel), and AGN_L100N512 (right 
panel), all at z = 0. Model WDENS uses an implementation of winds driven by supernovae that is more efficient in massive haloes than 
the prescription used in model REF. Model AGN uses the same recipe for supernova feedback as model REF, but also includes feedback 
from AGN. While there are some small differences, the trends are very similar for the three models. 



at successively higher redshiftqj, until we have exhausted 
all of them. For every particle we thus have values for the 
temperature, density, and metal mass as well as Zhaio and 
rahaio (unless no corresponding halo is found in our outputs). 
Note that the temperature and density recorded are for z — 
(or z = 2) , whereas rahaio is the mass of the host halo at 
redshift z h& i Q . 

If metals in the IGM were ejected from the haloes 
in which they formed shortly after they were released by 
stars, then the enrichment redshift, zz, and Zhaio should 
roughly agree. We would expect the two to be similar, as 
both metal enrichment and outflows are driven by massive 
stars, with Zhaio somewhat smaller than zz, simply because 
it takes some time for gas that has been ejected to reach the 
virial radius and leave the halo. However, even if enrichment 
and ejection happen nearly simultaneously, we would ex- 
pect differences of the order of the time scales corresponding 
to the time intervals between snapshots. We have checked 
this and found that Zhaio and z z do indeed agree very well 
for almost all temperatures and densities. This is true for 
both z = and z = 2. An exception is hot, dense gas 
(10 6 K <T< 10 7 K,p - 10 2 (p)) at z = 0, for which z haYo is 
significantly smaller than zz> This may be because this gas 
recently moved from the ICM across (our definition of) the 
halo boundary, but was enriched at relatively high redshift. 

The colour scale in Fig. [6] indicates the metal mass 
weighted mean log 10 rahaio as a function of the gas tem- 
perature and density for L050N512 at z = (left) and for 
L025N512 at z — 2 (right). Particles that were never part of 
a halo containing at least 100 dark matter particles were ex- 
cluded. These particles account for about 4 and 0.3 percent 
of the total metal mass in gas at z — and 2, respectively. 
Recall that gas with p > 10 2 (p) typically resides in a halo. 

There is a strong gradient with density in both panels. 
The lower the density, the lower the mass of the objects that 
were responsible for the enrichment. The enrichment of the 
low-density IGM (p < 10 (p)) is dominated by the smallest, 

5 Halo catalogues were generated at all simulation output red- 
shifts. These are spaced at Az = 0.125 for < z < 0.5, at 
Az = 0.25 for 0.5 < z < 4, and at Az = 0.5 for 4 < z < 9. 
We also have outputs for z = 10, 12.5, 15, and 20. 



resolved haloes, particularly at z = 0. This explains the rela- 
tively poor convergence that we found for these components 
in the previous section. Comparing with Fig. [2j we see that 
the low-density photo-ionised IGM (p < 10 (p), T < 10 45 K) 
is nearly completely missing from Fig. [6] This gas phase 
has therefore only been enriched by haloes below our res- 
olution limit and/or by intergalactic stars. The latter can 
occur when stars are ejected from galaxies via some dynam- 
ical process and then enrich the ambient IGM. 

There is a positive gradient of rahaio with T for high 
densities (p>10 2 (p)). Since most of these particles reside 
in haloes, this merely reflects the fact that the virial tem- 
perature is an increasing function of mass. 

Fig. [7] compares the reference model to model WDENS, 
which uses an implementation of galactic winds that makes 
them more efficient in high-mass haloes, and to model AGN, 
which includes AGN feedback. While there are some small 
differences, the trends are identical. Thus, the results are ro- 
bust to variations in the numerical implementation of galac- 
tic winds and to the inclusion of AGN feedback. Henceforth 
we will therefore only consider the reference model. 

Comparison of the left panel of Fig. [7] to the left panel 
of Fig. [6] which uses a mass (spatial) resolution that is eight 
(two) times higher, we see that the trends are also robust 
with respect to the numerical resolution. However, there are 
clear quantitative differences (note that the colour scales dif- 
fer), with the higher resolution simulation giving lower halo 
masses, particularly for masses near the resolution limit. 

Fig.[8]shows the metal mass weighted mean log 10 of the 
halo mass that last contained the metals as a function of den- 
sity (left panel) and temperature (right panel) for L050N512 
at z = (black) and for L025N512 at z = 2 (red). This plot 
therefore shows projections of Fig. [6] weighted by the metal 
mass in each pixel of that figure. The solid curves show the 
results for all the gas, while the other line styles correspond 
to cuts in the temperature (left) or density (right). Note 
that due to the finite resolution of the simulation, we may 
strongly overestimate the typical halo mass of origin when 
it is predicted to be of the order the minimum resolved halo 
mass. 

The solid curves in the left panel follow the dot-dashed 
curves at low densities (p<10 2 (p)) because most of the 
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Figure 8. Metal mass weighted mean log 10 rahalo? where rahalo is the total mass of the halo that last contained the metals, as a function 
of the gas density (left panel) or the gas temperature (right panel) at z = for model L050N512 (black) and at z = 2 for L025N512 
(red). The different line styles correspond to cuts in the temperature (left panel) or density (right panel) as indicated in the legends. 
The curves were cut when the number of particles per bin falls below 200. For both the diffuse IGM and the ISM, the typical halo mass 
that last contained the metals increases with the gas density. Gas at intermediate temperatures (T ~ 10 5 — 10 6 K) typically received its 
metals from lower mass haloes, but this merely reflects the fact that metals at these temperatures reside mostly in lower density gas 
than metals at lower or higher temperatures, as can be seen by comparing the different line styles in the right panel. 



1.0 
0.8 
0.6 

0.4 

0.2 
0.0 



ISM 





„ ^Diffuse IGM 
Cold Halo Gas 




9 10 11 12 13 149 

Log 10 (M hal0 [M ]) 



10 11 12 13 

Log 10 (M hal0 [M ]) 



14 9 10 11 12 13 

Log 10 (M hal0 [M ]) 



14 




10 11 12 13 14 9 

Log 10 (M hal0 [M ]) 



10 11 12 13 

Log 10 (M hal0 [M ]) 



10 11 12 13 

Log 10 (M hal0 [M ]) 



Figure 9. Cumulative probability distribution of the halo mass that last contained the metals, weighted by metal mass, at z = (top 
row) and z = 2 (bottom row) for stars and for the different gas phases defined in Fig. ^ The different colours correspond to different 
numerical resolutions. For metals that have never resided in a halo with more then 100 dark matter particles we set m^alo to the minimum 
resolved halo mass. The results are close to converged with respect to resolution except for the diffuse IGM and the WHIM. For these 
phases even our highest resolution runs may overestimate the halo mass. A wide range of halo masses contribute to the enrichment. 



metals residing in diffuse gas are hot (T > 10 4-5 K). Con- 
versely, for densities of order or greater than the star for- 
mation threshold, most of the metals are cold and the solid 
curves follow the dashed curves. We see similar behaviour 
in the right panel. 

The typical halo mass out of which the metals residing 



in the diffuse IGM were ejected increases strongly with the 
gas density. As the density increases above that typical of 
collapsed objects (p ~ 10 2 (p)), rahalo firsts declines, because 
cold metals begin to account for a substantial fraction of the 
metal mass and for lower mass haloes a greater fraction of 
the gas is in clumps. Beyond the star formation threshold 
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the mean halo mass increases again as more massive haloes 
sample higher pressure and thus, given our imposed equation 
of state for star- forming gas, higher density gas. The right 
panel shows that, for a fixed temperature, metals in higher 
density gas originate in higher mass haloes. The typical halo 
mass that last contained the metals reaches a minimum at 
intermediate temperatures (T ~ 10 5 — 10 6 K), because at 
these temperatures most of the metals reside at low-densities 
(p < 10 2 (p)), while high-density gas dominates the metal 
content at both low (T < 10 4 K) and high (T > 10 7 K) tem- 
peratures. 

In order to gain insight into the distribution of halo 
masses that contribute to the enrichment of a particular 
gas phase, we show the cumulative probability distribution 
function of rahaio in Fig. [9] Phrased differently, the figure 
shows the fraction of the metal mass that last resided in 
halo masses exceeding those plotted along the x-axis. The 
different columns and line styles correspond to the different 
regions in T — p space shown in Fig. [T] and the different 
colours show the results for different resolutions. The top 
and bottom rows are for z = and z = 2, respectively. 

For metals that never resided in a resolved halo we 
set rahaio to the dark matter mass of the smallest resolved 
haloes. The fraction of the metal mass ejected by unresolved 
haloes is greatest for the diffuse IGM at z = 0, for which it 
is about 12 percent. Comparing the two resolutions, we see 
that the results for the low-density phases, i.e. the diffuse 
IGM and the WHIM, are not converged. The situation is 
worst for the diffuse IGM at z = 0. Note that within these 
phases, convergence will be worse at lower densities as we 
have already demonstrated that the mean halo mass of ori- 
gin decreases with decreasing density. Fig. [9] clearly shows 
that a wide range of halo masses contributes to the enrich- 
ment of most gas phases. More than half of the metal mass 
in the diffuse IGM and the WHIM originated from haloes 
with total mass M < 10 11 M©, which corresponds to stellar 
masses smaller than 10 9 M . For the WHIM at z = 2 this 
fraction is about 80 %. As we noted in 33 if we separate 
the WHIM at p — 100 (p) (not shown), then the low-density 
component behaves similarly to the full WHIM, while the 
high-density component is more similar to the ICM. 



5 WHAT ENRICHED THE IGM AND WHEN? 

We conclude by simultaneously addressing the two questions 
that we investigated in the two previous sections: When was 
the IGM enriched and by what? Fig. [10] shows the metal 
mass distribution as a function of both zz and rahaio for 
model L050N512 at z — (left column) and for model 
L025N512 at z = 2 (right column). The top and bottom 
rows show the results for metals residing in the diffuse IGM 
and the WHIM, respectively. Metals that are currently part 
of a halo were excluded (note that this corresponds to a 
large fraction of the metals for the case of the WHIM) . As 
discussed before, the finite resolution of our simulations will 
cause us to underestimate the importance of haloes with 
masses near and below the minimum halo mass that we con- 
sider. 

Focusing first on the diffuse IGM at z = (the top- 
left panel of Fig. [T0|) . we see that there is a strong anti- 
correlation between the enrichment redshift and the mass of 



the halo from which the metals were ejected. Early enrich- 
ment is more strongly dominated by lower mass galaxies. For 
metals residing in the diffuse IGM at z — 2 (top-right panel) 
this is also the case, but the enrichment seems to have been 
dominated by objects with rahaio ~ 10 11 M© (which corre- 
sponds to a stellar mass ra* ~ 10 9 M©), which is about two 
orders of magnitude greater than the minimum halo mass 
we consider. The smaller role for the lowest mass haloes at 
z — 2 is also apparent from a comparison of the solid curves 
in the middle panels of Fig. [9j It appears that the resolution 
of the L025N512 simulation is nearly sufficient to model the 
enrichment of the diffuse IGM at z — 2. 

We must keep in mind, however, that the diffuse IGM 
is not a single phase. We have already demonstrated that 
zz and rahaio are strongly increasing resp. decreasing func- 
tions of the gas density (see Figs. [4] and [8]). For example, 
for the diffuse IGM at z — 2 the metal mass weighted mean 
log 10 rahaio decreases from w 12 for an overdensity of 10 2 to 
~ 10 (i.e. too close to the resolution limit to be converged) 
at the mean density. 

It is likely that a significant fraction of the metal mass 
that resides at the high-density end of the diffuse IGM at 
z = 2, which is typically ejected by more massive haloes, 
will have fallen back into haloes by z — and will therefore 
no longer be part of the diffuse IGM. This effect can help 
explain why rahaio is smaller at z — than at z — 2 and why 
there is an ant i- correlation between zz and rahaio- Evolution 
of the halo mass function is likely also an important part of 
the explanation. At higher redshift s the masses of the haloes 
producing metals are typically lower. 

The enrichment history of the WHIM (bottom row of 
Fig. I10p shows some interesting differences from that of the 
diffuse IGM. At z = the ant i- correlation between zz and 
rahaio is much weaker. Relatively recent enrichment by low- 
mass haloes is more important for the WHIM. Presumably, 
metals that are ejected recently are shock-heated by the 
winds that carry them and many of them may not have 
had sufficient time to cool down to the temperatures of the 
diffuse IGM. High-mass haloes are also more important for 
the WHIM. This may largely be because high-mass haloes 
tend to be embedded in gas that is heated by gravitational 
accretion shocks and is therefore sufficiently hot to be con- 
sidered part of the WHIM. 

The results for the other OWLS models considered here 
are very similar (not shown), although higher mass haloes 
are slightly more important for WDENS and AGN, presum- 
ably because in these models feedback is more efficient in 
such objects. 



6 CONCLUSIONS 

We have used cosmological, hydrodynamical simulations 
that include radiative cooling, star formation, stellar evo- 
lution and supernova feedback, and, in some cases, AGN 
feedback, to investigate the metal enrichment history of the 
Universe. Specifically, we took advantage of the Lagrangian 
nature of our simulations to investigate when gas that ends 
up in a particular gas phase at redshift z = or z = 2 re- 
ceived its metals and what the masses were of the haloes 
that last contained these metals. 

We considered different physical models, all taken from 
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Figure 10. Distribution of metal mass as a function of ra^alo? the metal mass weighted mean mass of the last resolved halo that contained 
the metals, and zz, the metal mass weighted mean enrichment redshift. Results are shown for the diffuse IGM (top row) and the WHIM 
(bottom row) at z = for L050N512 (left column) and at z = 2 for L025N512 (right column). The contours are logarithmically spaced 
by 0.5 dex. Metals that are currently in a halo and metals that have never been part of a resolved halo were both excluded. Older metals 
were typically ejected by lower mass haloes, particularly if they end up in the diffuse IGM. The WHIM was enriched by a wide range of 
halo masses. 



the OWLS project, including simulations with a different 
implementation of supernova feedback and a model that in- 
cludes AGN feedback. While there are some small differ- 
ences, such as slightly higher enrichment redshifts for models 
with more efficient feedback in massive galaxies, the main 
trends are strikingly similar for all models. We cannot rule 
out that some of our conclusions may depend on our pre- 
scriptions for feedback or that they are specific to SPH. 
However, the fact that three drastically different feedback 
models give such similar results suggests that our conclu- 
sions are robust. 

The time since the enrichment varies most strongly with 
the gas density. Metals in lower density gas are typically 
much older than metals in high-density gas, a trend that 
extends over ten orders of magnitude in gas density. At least 
half of the metals that reside in the diffuse IGM at z — 
(z — 2) were ejected above redshift two (three). 

The enrichment redshift also varies with the tempera- 
ture of the gas, but this mostly reflects the fact that the 
typical gas density of enriched gas varies with its tempera- 
ture. For gas that is dense enough to be able to cool, the time 
since enrichment correlates well with the cooling time. Gas 
with shorter cooling times received its metals more recently. 



The typical mass of the haloes from which the metals 
residing in the IGM were ejected, increases rapidly with the 
gas density. At least half of the metal mass was ejected by 
haloes with total masses less than 1O 11 M0, which corre- 
sponds to stellar masses smaller than 10 9 M©. For the low- 
redshift IGM the mass of the dominant IGM polluters may 
be substantially smaller than this because our predictions 
for present day low- density gas are not close to converged 
with respect to the numerical resolution. 

The age of the metals ending up in the diffuse IGM is 
strongly anti-correlated with the mass of the haloes from 
which they were ejected. In other words, older metals were 
typically ejected by lower mass haloes. This ant i- correlation 
is less strong for the hotter part of the IGM (the WHIM), 
for which recent enrichment by low-mass haloes is more im- 
portant. While massive haloes (> 10 12 M ) are unimportant 
for the enrichment of the diffuse IGM, they do contribute to 
the pollution of the warm- hot IGM at low redshift. 

Our main and strongest conclusion is that metals resid- 
ing in lower density gas were typically ejected earlier and 
by lower mass haloes. This suggests that travel-time is a 
li miting factor for the e nrichment of the IGM, as proposed 
bv lAguirre et al] (|2001br ). Metals that have been ejected re- 
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cently simply have not had sufficient time to reach the low- 
density gas far from galaxies. The ant i- correlation between 
metal age and halo mass would then follow because typical 
galaxy masses are lower at higher redshifts. Another mecha- 
nism that could explain the trends is fall back onto galaxies. 
If metals ejec ted by lower mass galaxi es are less likely to 
fall back (e.g. lOppenheimer et al.l 120091 ) or if they typically 
leave the haloes with larger velocities, then we would ex- 
pect lower mass galaxies to be more important for gas that 
is further away. Finally, low-mass galaxies are more weakly 
clustered than high-mass galaxies, thus low-density gas is 
more accessible to low-mass galaxies. 

It is not immediately obvious how one could compare 
this result with observations. One way would be to consider 
abundance ratios of elements in various phases as a func- 
tion of density and redshift in an attempt to find a signa- 
ture that could be searched for in high-redshift observations. 
However, uncertainties in the nucleosynthetic yields and the 
type la supernova rates (see Wiers ma et al.ll2009br L as well 
as in the ionization corrections, may make such a compar- 
ison troublesome. Correlations between absorbers and the 
distances to different types of galaxies may also help shed 
light on the en richment mechan ism (e.g . [Steidel et al. 2010; 
Che n fe Mulchaevl 120091 : IWilman et al.l 12007: St ocke et all 



2006), although is is important to keep in mind that the 



metals that are observed to be near a gal axy may have been 
injec t ed by its lower-ma ss progenitors (|Porciani fc Madaul 
120051 ; ISca nnapiecol 120051 ). In fact, our results would imply 
that this is likely to be the case for metals in low-density 
gas. We leave an analysis of abundance ratios and other ob- 
servational comparisons to a future work. 

The importance of low-mass galaxies is a challenge for 
simulations of the enrichment of the IGM, as it implies the 
need for high-resolution in large volumes. On the other hand, 
our results suggest that the low-density IGM provides us 
with an exciting opportunity to study the consequences of 
feedback in low-mass galaxies, a key but poorly understood 
ingredient of models of galaxy formation, and that it pro- 
vides us with a fossil record of galaxy formation in the high- 
redshift Universe. 
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